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We consider an approximation procedure to evaluate the finite-temperature one-loop fermionic 
density in the presence of a chiral background field which systematically incorporates effects from 
inhomogeneities in the chiral field through a derivative expansion. We apply the method to the 
case of a simple low-energy effective chiral model which is commonly used in the study of the chiral 
phase transition, the linear (T-model coupled to quarks. The modifications in the effective potential 
and their consequences for the bubble nucleation process are discussed. 



I. INTRODUCTION 



o 
o 

CsJ ' It is commonly accepted that QCD at sufficiently high temperatures undergoes a phase transition to a new state of 

matter, the quark-gluon plasma (QGP), which was presumably present in the early universe [l|,l3l- CompeUing lattice 

^^ QCD results corroborate this belief .3], and experiments in ultra-relativistic heavy-ion collisions |j,|a] at BNL-RHIC 

V^ , have recently shown data that clearly point to a new state of matter |^ . 

r — ' To model the mechanism of chiral symmetry breaking present in QCD, and to study the dynamics of pha se con- 

,—1 I version after a temperature-driven chiral transition, one can resort to low-energy effective models 0, H S Il3 IHl u2, 
Ha. I14L I15L llq. In particular, to study the mechanisms of bubble nucleation and spinodal decomposition in a hot 

CN ' expanding plasma [13, it is common to adopt the linear cr- model coupled to quarks ISj, where the latter comprise 
K*" . the hydrodynamic degrees of freedom of the system. The gas of quarks provides a thermal bath in which the long- 
wavelength modes of the chiral field evolve, and the latter plays the role of an order parameter in a Landau-Ginzburg 
approach to the description of the chiral phase transition 11, 12, 13, 14]. The standard procedure is then integrating 

^P over the fermionic degrees of freedom, using a classical approximation for the chiral field, to obtain a formal expression 

f^ , for the thermodynamic potential: 



5' r!(r,^,0) = F(0)-^lndet 



(G^i + M(0)) 



T 



(1) 



II ' where V{(f>) is the classical self-interaction potential for the bosonic sector, Ge is the fermionic Euclidean propagator, 
Oni M((j)) is the effective fermion mass in the presence of the chiral field background, T is the temperature and V is the 
2^ ■ volume of the system. From the thermodynamic potential (^), one can obtain all the physical quantities of interest. 
'■^ ' To actually compute correlation functions and thermodynamic quantities, one has to evaluate the fermionic de- 

terminant that results from the functional integration over the quark fields within some approximation scheme. 
Alternatively, one can consider the fermionic density which will appear as a source term in the equation of motion 
for the chiral field. In the case of one-dimensional systems one can often resort to exact analytical methods, such as 
C^ ' the inverse scattering technique |19| |. In practice, however, the determinant is usually calculated to one-loop order 
assuming a homogeneous and static background field 20;j. Nevertheless, for a system that is in the process of phase 
conversion after a chiral transition, one expects inhomogeneities in the chiral field configuration due to fluctuations to 
play a major role in driving the system to the true ground state. Hence, their effects should in principle be included 
in the computation of the fermionic determinant. 

In the case of high-energy heavy ion collisions, hydrodynamical studies have shown that significant density inho- 
mogeneities may develop dynamically when the chiral transition to the broken symmetry phase takes place | l3| (see 
also [IjI for an analysis in a different context). Their pattern and intensity might indeed provide some insight on 
the nature of the transition as well as on the location of an eventual critical point. If the freeze-out in heavy ion 
collisions occurs shortly after a first-order chiral transition, inhomogeneities generated during the late stages of the 
nonequilibrium evolution of the order parameter might leave imprints on the final spatial distributions and even on 
the integrated, inclusive abundances I21j . 

In this paper we consider an approximation procedure to evaluate the finite-temperature fermionic density in the 
presence of a chiral background field which systematically incorporates effects from inhomogeneities in the bosonic 
field through a gradient expansion. The method is valid for the case in which the chiral field varies smoothly, and 
allows one to extract information from its long-wavelength behavior, incorporating corrections order by order in the 
derivatives of the field 22] . This approach has been successfully used to treat systems of low-dimensionality at zero 
temperature in condensed matter physics 23]. Here we consider a three-dimensional system at finite temperature. We 
apply the method to the case of the linear cr-model coupled to quarks, which provides a convenient framework for the 



study of bubble nucleation and spinodal decomposition in the case of a first-order chiral transition. Nevertheless, the 
results presented below are quite general and may be of interest also in cosmology or in condensed matter systems. 

The paper is organized as follows. Section II presents briefly the low-energy effective model adopted in this paper. 
In Section III we introduce the method to incorporate systematically effects from inhomogeneities in the chiral field 
in the computation of the fcrmionic density. Results for the (well-known) leading term and for the first non-trivial 
corrections are discussed in Section IV. There, we also consider the modifications undergone by the effective potential 
and their consequences to the process of nucleation. Section V contains our final remarks. 

II. EFFECTIVE MODEL 

Let us consider a scalar field (f) coupled to fermions V' according to the Lagrangian 

£ = i^li^d^ + M7" - M{^)]^ + ^d^<l> d^'cj) - y(</>) . (2) 

where /i is the fermionic chemical potential, M(0) is the effective mass of the fermions and V{(j)) is a self- interaction 
potential for the bosonic field. 

In the case of the linear cr-model coupled to quarks, (f> represents the a direction of the chiral field <& — (cr, tt), where 
TT* are pseudoscalar fields playing the role of the pions, which we drop here for simplicit y. T he pion directions play no 
major role in the process of phase conversion we have in mind, as was argued in Ref. [ll|, so we focus on the sigma 
direction in what follows. However, the coupling of pions to the quark fields might be quantitatively important in the 
computation of the fermionic determinant inhomogeneity corrections. This issue makes the computation technically 
more involved and will be addressed in a future publication. The field ip plays the role of the constituent-quark 
field q — {u,d), and ^ = /J,g is the quark chemical potential. The "effective mass" is given by M{(j)) — g\a\, and 
V{^) = (A^/4)(f7^ -I- 7?^ — v'^Y ~ ^q'^ is the self-interaction potential for <I>. The parameters above are chosen such 
that chiral SUl{2) (8) SUr{2) symmetry is spontaneously broken in the vacuum. The vacuum expectation values of 
the condensates are (cr) = /jr and (tt) = 0, where /^r = 93 MeV is the pion decay constant. The explicit symmetry 
breaking term is due to the finite current-quark masses and is determined by the PCAC relation, giving hq = /jrmj, 
where m^r = 138 MeV is the pion mass. This yields v^ — f"^ — rn^/X^. The value of A^ = 20 leads to a cr-mass, 
TO^ = 2A^/^ -I- m^, equal to 600 MeV. In mean field theory, the purely bosonic part of this Lagrangian exhibits a 
second-order phase transition J2J| at Tc — \f2v if the explicit symmetry breaking term, hq, is dropped. For hq ^ 0, the 
transition becomes a smooth crossover from the restored to broken symmetry phases. For g > 0, one has to include a 
finite-temperature one-loop contribution from the quark fcrmionic determinant to the effective potential as indicated 
in Eq. Q. When the coupling between quarks and the chiral field, g, is large enough, the system exhibits a first-order 
phase transition even at /z = [allBliSl- When we decrease g, the strength of this first-order transition is weakened. 
At g ~ 3.7, the latent heat vanishes and we have a second-order critical point at /U = 0. In what follows we keep 
the explicit symmetry breaking term hqa and consider the case g = 5.5, where the first-order line goes all the way 
down to /i = 0, since we are mainly concerned with the effects from inhomogeneities in the process of homogeneous 
nucleation. 

The Euler-Lagrange equation for static chiral field configurations contains a term which represents the fcrmionic 
density, p : 

dV 

'^^(t^ = ^+gp{T,fi,^), (3) 



and the density of fermions at a given point Xq has the form 



p{xo) = Sp(xo 



G-^^ + M{x) 



xo) , (4) 



where \xo) is a position eigenstate with eigenvalue xq, and Sp represents a trace over fermionic degrees of freedom, 
such as color, spin and isospin. 

Assuming a homogeneous background field, one can compute the one-loop fermionic density in a simple way |2(]| |. In 
this case, the correction coming from the integration over the fermions can be directly incorporated into an effective 
potential for the chiral field, as will be shown below. However, perfect homogeneity is a very strong hypothesis if 
one is interested in the dynamics of a phase transition. On the other hand, the correct determinant would have to 
be computed with an arbitrary profile for the background field. In a few examples, one can do it formally for one- 
dimensional systems 19]. For higher dimensions, however, one must adopt some approximation scheme to take into 



account inhomogeneity effects. In tfie next section, we present a framework to incorporate systematically derivative 
corrections to the density p{x). The only assumption made on the behavior of the background field is that it varies 
very smoothly. 



III. INHOMOGENEITY CORRECTIONS 



In order to take into account inhomogeneity effects of the chiral background field, 0, 
dependence of M in Q , we resort to a derivative expansion as explained below. 

In momentum representation, the expression for the fermionic density assumes the form 



p{xo) ^ SpT^ 



<fik 

(27r)3 



-ik-Xf) 



1 



Ak-XQ 



-1^{iLUn + /Lt) - 7 • fc + M{x) 



encoded in the position 



(5) 



where ojn = (2n + 1)'kT are Matsubara frequencies for fermions [20j. One can transfer the xq dependence to M{x) 
through a unitary transformation, obtaining 



p(^o) 



''^Vlwi'Vo 



7"(ia;„ + ^) - 7 • fc + M{x + xq) 



(6) 



where one should notice that xq is a c-number, not an operator. 
Now we expand M{x + xq) around xq: 

M(X + fo) = M{xq) + AM(X, Xq) = 



M{xo) + ViM{xo)x' + -\7i\7jM{xo)x'x^ + 



(7) 



and use a;* 



-iVz- to write 



piS,) = SpTY^J^- 



(27r)3 -^o(j^^^ + ^) _ ^ . /j + M{xo) 



l + AM{-iVk,,xo] 



-f^iiujn + /i) - 7 • fc + M{xo) 



(8) 



To study the dynamics of phase conversion after a chiral transition, one can focus on the long-wavelength properties 
of the chiral field. From now on we assume that the static background, M{x), varies smoothly and fermions transfer 
a small ammount of momentum to the chiral field, so that AM/M << 1. Under this assumption, we can expand the 
expression inside brackets in Eq. (|HJ| in a power series: 



pi:g)=SpTj2j^- 



1 



(27r)3 ^0(j^^^ ^^) _ ^ . fc _,_ ^/(^) 
Eq. ® , together with 



E(-i)^ 



AM(-iVfc,,f) 



-f°{iujn + fJ-)-^-k + M{x) 



(9) 



AM{~iVki,x) = VjM(x) 



'k. 



-V,VjM(f) 



Vfc^Vfc^. 



(10) 



provides a systematic procedure to incorporate corrections brought about by inhomogeneities in the chiral field to the 
quark density, so that one can calculate p{x) — po{x) + pi{x) + P2{x) + • • • order by order in powers of the derivative 
of the background, M{x). 

The new corrections will bring higher-order derivatives to the equation of motion for the chiral field. In particular, 
as will be seen below, the first non-trivial inhomogeneity contribution will modify the Laplacian term in Eq. ^, and 
can be seen as a correction to the surface tension in the process of bubble nucleation. 

This is a quite general method to approximate the fermionic density and could be used in a variety of low-energy 
effective field theory models for the study of the dynamics of the chiral transition. In the next section we apply this 
method to the case of the linear cr-model coupled to quarks. 



IV. RESULTS 



A. Leading term 



The leading-order term in this gradient expansion for p{x) can be calculated in the standard fashion i2C)l| and yields 
the well-known mean field result for the scalar quark density 



Po 



(27r)3 e[EU<t>)-f^,]/T + 1 



+ {Pq -^ -fJ-q) , 



(11) 



where Vq = 12 is the color-spin-isospin degeneracy factor, Ek{4>) = (P -I- M^{(j))y/^, and M{(j)) = g\(j)\ plays the role 
of an effective mass for the quarks. The net effect of this leading term is correcting the potential for the chiral field, 
so that we can rewrite Eq. Q as 



V"(i) = 



dVe 



eff 



(12) 



where Veff = V{(j>) + Vq{(j>) and 



V„ 



-VqT 



ffk 



(2 



^ln(e[^'=(*)-'''l/^ + l)+(A.,--M,) 



(13) 



The potentials Vq and Veff for several values of the temperature (at ^iq = 0) are displayed in Fig. 1 and Fig. 2, 
respectively, assuming 5 = 5.5. In this case, the chiral phase transition is of first order even for a vanishing chemical 
potential, an d Tr « 124 MeV. The barrier for nucleation disappears at T^p k. 108 MeV, where the system reaches the 
spinodal line [ij. 
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FIG. 1: Vg(<^) for different values of the temperature at /iq = and for g = 5.5. 



This kind of effective potential is commonly used as the coarse-grained thermodynamic potential in a phenomeno- 
logical description of the chiral transition for an expanding quark-gluon plasma created in a high-energy heavy-ion 
collision (ol. HoI HiL ri3| . From the modified field equation (|12ll . one can study, for instance, the phenomena of bubble 
nucleation and spinodal decomposition. However, the presence of a non-trivial background field configuration, e.g. a 
bubble, can in principle dramatically modify the Dirac spectrum |25j, hence the determinant. In the case of condensed 
matter systems, where electronic doping often plays a major role, the presence of fermionic bound states can deeply 
affect the dynamics of the phase transition. This is the case in the presence of a bubble background, where besides 
unstable critical bubbles one can find metastable configurations, dep ending on the relative occupation of bound states, 
and a modification in the value of the nucleation rate [l9l l2a. l27j . In the case of a chiral model, analogous effects 
can in principle appear for a nonzero quark chemical potential. In any case, one expects the effective potential to be 
modified by the effect of fluctuations of the chiral field on the fermionic density, which motivates the investigation of 
the next term in the expansion, which contains some information about the inhomogeneity of the bosonic field. 
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FIG. 2: Veff{(fi) in the a direction for different values of tfie temperature at Hq = and for g — 5.5. 



B. First corrections 

The next non-trivial term in the expansion contains two contributions: one coming from pi and another from p2- 
This is due to the rearrangement of powers of the gradient operator. This term will correct the Laplacian piece in the 
chiral field equation. Dropping zero-temperature contributions which can be absorbed by a redefinition of the bare 
parameters in V, a long but straightforward calculation yields 



{Pl + P2) = -i^^^M) Wg{T, fig,Cl)) , 



(14) 



where 
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(16) 



and np{Ek, p,q) is the Fermi-Dirac distribution. The derivation of H{Ek, T, p,q) is not particularly illuminating (a few 
steps are presented in the appendix). However, in the low-temperature limit, corresponding to /3M >> 1, the integral 
above is strongly suppressed for high values of k, and the leading term has the much simpler form 



Wq{T,Pq,^)^^^e-^\''^\m\?" ^ 



(17) 



which gives a a better idea of the profile of the first inhomogeneity correction. One can already anticipate that it will 
be concentrated in the same region where the homogeneous correction was significant, i.e. /3M < 1 (cf. Fig. 1), being 



T=50 MeV 

T= 100 MeV 

T=150MeV 




150 

^ (MeV) 



FIG. 3: gWq{(j>) for different values of the temperature at ^g = and for g = 5.5. 



exponentially suppressed for higher values of the field. In fact, a numerical study of the complete Wq shows that this 
function is peaked around = and non-negligible for /3|(^| < 1 (see Fig. 3). 

The Euler-Lagrange equation for the chiral field up to this order in the gradient expansion reads 



V^^^[l+gWgiT,^lg,^)]- 



dVeff[T,llq, 

d(j) 



VIM) 



(18) 



Here we used the fact that (l+gWg) is a positive definite quantity, and defined a new "effective potential" that contains 
all the corrections up to this order in the gradient expansion, Vm- One should not confuse Vjn with the standard 
definition of the one- loop effective potential in field theory derived for a constant background ^8|. Nevertheless, we 
keep the name effective potential for Vm for convenience in the description of nucleation that follows. 

The complete new effective potential can be obtained from our previous results by numerical integration. In order to 
proceed analytically, though, we choose to fit its derivative, which we know exactly up to this order, by a polynomial 
of the fifth degree. Actually, we know that the commonly used effective potential, T4// = V +J4, can hardly be 
distinguished from a fit with a polynomial of sixth degree in the region of interest for nucleation '2!|| . Working with 
fits will be most convenient for using well-known results in the thin- wall approximation to estimate physical quantities 
that are relevant for nucleation, such as the surface tension and the free energy of the critical bubble. Results for the 
fits of V'/„(0) for fiq = are shown in Fig. 4. 
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FIG. 4: Fits for Vi'„{<j>) for different values of the temperature at fiq = and for g — 5.5. The lines shown in the box correspond 
to the exact results. 



We can now integrate analytically the polynomial approximation to the derivative of the complete effective potential. 
In Fig. 5 we display the curves for Ve/f and Vm for a few values of temperature and Hq = 0. 
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FIG. 5: Veff and Vin{(l>) for different values of the temperature T = [108 (a), 116 (b), 124 (c), 132 (d)] MeV at ^^ = and for 
<7 = 5.5. 

From Fig. 5 one can notice a few consequences of the inhomogeneity correction. The first general effect is the 
smoothening of the effective potential. In particular, and most importantly, the barrier between the symmetric phase 
and the broken phase is significantly diminished, as well as the depth of the broken phase minimum, although we still 
have a first-order phase transition barrier. Therefore, one can expect an augmentation in the bubble nucleation rate. 
In principle, one should have better results from calculations within the thin-wall approximation. Also, the critical 
temperature moves up slightly. 

C. Effects on nucleation 

Let us now consider the effects of the first inhomogeneity correction on the process of phase conversion driven by 
the nucleation of bubbles [30|]. To work with approximate analytic formulas, we follow Ref. |ll| and express Vin over 
the range < < T in the familiar Landau-Ginzburg form 



Vr 



off 



J2 «" *' 



(19) 



Tl = 



Although this approximation is obviously incapable of reproducing all three minima of Vin , this polynomial form is 
found to provide a good quantitative description of Vm in the region of interest for nucleation, i.e. where the minima 
for the symmetric and broken phases, as well as the barrier between them, are located. 



A quartic potential such as Eq. (|19|l can always be rewritten in the form 

U{v)^a{>/-ay+j^ . (20) 

The coefficients above are defined as follows: 

(21) 

(22) 

(23) 

(24) 

The new potential U{ip) reproduces the original Vi„((^) up to a shift in the zero of energy. We are interested in the 
effective potential only between Tc and Tsp. At T^, we will have two distinct minima of equal depth. This clearly 
corresponds to the choice j = in Eq. H20() so that U has minima at (^ = ±a and a maximum ai ip — 0. The minimum 
at Lp = —a and the maximum move closer together as the temperature is lowered and merge at Tgp. Thus, the spinodal 
requires j/ao? — — 8/3\/3 in Eq. if^ . The parameter jjac? falls roughly linearly from 0, at T = T^ to — 8/3\/3 at 
the spinodal. 

The explicit form of the critical bubble in the thin-wall limit is then given by |l9| 
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where iff is the new false vacuum, Re is the radius of the critical bubble, and ^ 
a measure of the wall thickness. The thin-wall limit corresponds to £,/Rc <C 1 



(25) 

2/to, with m^ = U"{(pf), is 
I, which can be rewritten as 



(3|j|/8q;o'^) <C 1. This small parameter has the value of l/V^ at the spinodal, which suggests that the thin-wall 
approximation might be qualitatively reliable for our purposes. Nevertheless, it was shown in |ll| that the thin- 
wall limit becomes very imprecise as one approaches the spinodal. In this vein, the analysis presented below is to 
be regarded as semi-quantitative. To be consistent we compare results from the homogeneous calculation to those 
including the inhomogeneity correction within the same approximation. 
In terms of the parameters a, o, and j defined above, we find 
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in the thin- wall limit. Determination of the critical radius requires the surface tension, S, defined as 

fdipb\\^ 2 



dr , , 

V dr / 



3a^-^ 



(26) 
(27) 

(28) 



The critical radius then becomes Re = (2S/AC/), where AC/ = V{(j)f) — V{(j)t) « 2o|j|. The free energy of a critical 
bubble is finally given by Fi, — (47rl]/3)i?^. From knowledge of Fb, one can evaluate the nucleation rate T ^ e^^^l^ . 
In calculating thin- wall properties, we shall use the approximate forms for (/jj, (/jy, S, and At/ for all values of the 
potential parameters. 

To illustrate the effect from the inhomogeneity correction, we compute the critical radius and F^/T for three 
different values of the temperature. For T = 108 MeV, corresponding to the spinodal temperature, which is not 
modified by the first inhomogeneity correction, the corrected values are Re ~ 0.98 fm and F^/T « 0.20, as compared 
to Re ~ 1.1 fm and F^/T « 0.9 in the homogeneous case. The same computation for T — 116 MeV yields Re ~ 2.15 fm 
and Fb/T « 1.14, as compared to Re « 2.2 fm and Fb/T « 2.1. At T = 124 MeV, which corresponds to the critical 
temperature for the homogeneous case, the critical radius and Fb/T diverge in the homogeneous computation, whereas 
Re ~ 35 fm and Fb/T « 394 including inhomogeneities. The numbers above clearly indicate that the formation of 
critical bubbles is much less suppressed in the scenario with inhomogeneities, which will in principle accelerate the 
phase conversion process after the chiral transition. 



SUMMARY AND OUTLOOK 



We have introduced a systematic procedure to evaluate inhomogeneity corrections to the finite-temperature 
fermionic density in the presence of a chiral background field, which incorporates effects from fluctuations in the 
bosonic field through a gradient expansion at finite temperature and density. Higher-order contributions give more 
non-local corrections to the effective Euler-Lagrange equation for the chiral field, and the condition for the validity 
of the method is a smooth variation of the chiral field, which should be enough in the analysis of its long-wavelength 
behavior in the phase transition. 

Incorporating the first inhomogeneity correction in the computation of the effective potential of the linear cr-model 
coupled to quarks, we found that the latter is significantly modified. Besides a general smoothening of the potential, the 
critical temperature moves upward and the hight of the barrier separating the symmetric and the broken phase vacua 
diminishes appreciably. As a direct consequence, the radius of the critical bubble goes down, as well as its free energy, 
and the process of nucleation is facilitated. Although the numbers presented above should be regarded as simple 
estimates, since they rely on a number of approximations, the qualitative behavior is clear. In a detailed quantitative 
analysis, one should not only integrate numerically the effective potential and relax the thin-wall approximation, but 
also include the pion-quark interaction in the computation of the fermionic density. We believe that the contribution 
from the pion sector will enhance the effect from inhomogeneities. 

In all the discussion above, we intentionally ignored corrections coming from bosonic fluctuations, which would 
result in a bosonic determinant correction to the effective potential j^ll] . To focus on the effect of an inhomogeneous 
background field on the fermionic density, we treated the scalar field essentially as a "heavy" (classical) field, whereas 
fermions were assumed to be "light" . 

Experimental signatures of inhomogeneities for high-energy heavy ion collisions were discussed, for instance, in 
Refs. [i3i|2ll- In particular, inhomogeneities seem to favor an "explosive" scenario jig for the phase conversion even 
at early stages of nucleation. However, one should first incorporate dissipation and noise effects, which tend to retard 
the explosion ,29.] . before estimating the time scales involved. This analysis will be left for a future publication. 



Acknowledgments 

The authors are grateful to A. Dumitru for a critical reading of the manuscript and several suggestions. We also 
thank D.G. Barci, H. Boschi-Filho, C.A.A. de Carvalho and T. Kodama for discussions. This work was partially 
supported by CAPES, CNPq, PAPER J and FUJB/UFRJ. 

APPENDIX A 

In this appendix we sketch the main steps to build the function Wq{T, fj,q, (f>) that corrects the Laplacian in the 
Euler-Lagrange equation for the chiral field. 

The first inhomogeneity correction has contributions from pi and p2. The contribution coming from pi is propor- 
tional to ViV jcj): 



,.,., . _.., /^ I^^St^ (iv,V,M<..,..) P^f^^ 



where we use a compact notation for the sum-integrals 



K 



(2^ 



(Al) 



(A2) 



K — (fco, k) = (iujn — fJ-qTk), and Tr^ is a trace over Dirac gamma matrices. 
There is also a contribution proportional to \Ji(j)V j(j) coming from p2'- 



Pi\x) ^Tr^ 



7°fco -7- k- M{x) 



K 



kl - El 



Up to this order in derivatives, we can write 
d?k 



ViM(f)i' 



7°fco - 7 • fc - M{x) 



kl - El 



(A3) 



Pi 



/d k r - -> ^ 

-^^Vq[^HI{x)[F{Ek,^iq,P)+F{EJ,,~^Xq,P)]^{k■^fM{x)[G{Ek,^Jiq,P)+G{Ek,-pq,P)]\ 
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(A4) 



where /3 = 1/T, 

-/3 



and 



1 M2\ /e2'3(^k-A'9) -e'^(^'="^<') 
"24"e|"J (^ (e/J(i^fc-A',) + 1)3 

1 3 M2\ / e'^(^'=-''') 



1 3 M2\ / 1 



23^3 24^5.^ \^(g/3(i3.-p,) + l) 



(A5) 



2'^ El ) \ (g/J(£fc-M,) + 1)5 



3 / 1 M^ 5M4\ /e^'^^^""^") -4e2'^(^'«"^9) +e'^(^'«"^9) 



23 £^4. 2E%) \ (e/5(^'=-A',) + 1)4 

1 3 M^ 45M4\ /e2'3(-^'=-Ai9) _ e/3(Bfc-M,) 



'^ ly 23£;3. ' 22£;5. ' 2'^E]^) \ (e/3(Bfc-/^,) + 1)3 

^ 3 / 1 5M2 M4\ / e'3(s»=-A*,) 

-/37T7 -:T7^ + 7T7Fr+35- 



22 V 2£4_ ' 2E% E%) \{ePiEk-i'o) + lY 

1 5M2 M4\ / 1 \ 

35-^ 7-^7^-7-^^ . (A6) 



22 V 2E% 2EI. E^j^J yie'^iEj.-f^,) + 1) 

Rewriting the term proportional to fc • V in a more convenient form, using i?| ^ k'^ + M'^ , and exploring some 
symmetries in the integrands, it is straightforward to arrive at the final form: 

Pi+P2 = v20(f)0y"dfcP[i7(Sfc,Ai„/3)+if(£fc,-A*9,/?)] , (A7) 



where 



H{Ek,^ig,P) = F{Ek,^ig,P) + h''G{Ek,^ig,l3) , (A8) 



which gives the correction to the Laplacian. 
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